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ABSTRACT 

We have studied the velocity field of the Blue Compact Dwarf galaxy Mrk86 
(NGC 2537) using data provided by 14 long-slit optical spectra obtained in 10 dif- 
ferent orientations and positions. This kinematical information is complemented with 
narrow- band ([OIII]5007A and Ha) and broad-band (B, V, Gunn-r and K) imaging. 
The analysis of the galaxy global velocity field suggests that the ionized gas could 
be distributed in a rotating inclined disk, with projected central angular velocity of 
SI = 34 kms _1 kpc _1 . The comparison between the stellar, HI and modeled dark 
matter density profile, indicates that the total mass within its optical radius is dom- 
inated by the stellar component. Peculiarities observed in its velocity field can be 
explained by irregularities in the ionized gas distribution or local motions induced by 
star formation. 

Kinematical evidences for two expanding bubbles, Mrk86-B and Mrk86-C, are 
given. They show expanding velocities of 34kms _1 and 17kms _1 , Ha luminosities of 
3xl0 38 ergs _1 and 1.7xl0 39 ergs _1 , and physical radii of 374 and 120 pc, respectively. 
The change in the [SII]/Ha, [Nil] /Ha, [OII]/[OIII] and [OIII]/H/3 line ratios with the 
distance to the bubble precursor suggests a diminution in the ionization parameter and, 
in the case of Mrk86-B, an enhancement of the shock-excited gas emission. The optical- 
near-infrared colours of the bubble precursors are characteristic of low metallicity star 
forming regions (^0.2 Zq) with burst strengths of about 1 per cent in mass. 
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1 INTRODUCTION 

The study of the global velocity field of Blue Compact Dwarf 
galaxies (BCD hereafter; Thuan & Martin 1981) provides 
important clues about their gravitational potential, since 
these systems are rotationally supported (e.g. van Zee et 
al. 1998). 

High spatial resolution HI observations have shown 
that the rotation curves of Blue Compact Dwarf galaxies 
(Meurer, Staveley-Smith & Killeen 1998; van Zee et al. 1998) 
and dwarf irregulars (dl hereafter; Moore 1994; Flores & Pri- 
mack 1994) are nearly flat in the galaxy outer regions and 
have nearly constant velocity gradients within their optical 
radius. Also, optical studies of the velocity field of the ion- 
ized gas in BCDs obtain constant velocity gradients, char- 
acteristic of a solid-body rotation law (see Petrosian et al. 
1997 for I Zw 18). 

Although the neutral (see van Zee et al. 1998 and ref- 
erences therein) and molecular hydrogen (Young & Knezek 
1989; Israel, Tacconi & Bass 1998) are quite abundant in 
Blue Compact Dwarfs and dwarf irregulars, they are not 
enough to reproduce the flattening of the rotation curve. 



Like in spiral galaxies, the existence of this flattening in the 
rotation curve of dwarf galaxies has been related with the 
presence of large amounts of dark matter in galaxy outer 
regions (Carignan & Freeman 1988; Carignan & Beaulieu 
1989; Broeils 1992). The dark matter content derived indi- 
cates that dark matter is even more abundant in dwarfs that 
in more massive galaxies (see Moore 1994 and references 
therein). In fact, standard cold dark matter (CDM here- 
after) models predict that low-mass halos are denser than 
more massive systems, because their higher formation red- 
shift (Navarro, Frenk & White 1997, NFW hereafter). The 
density profiles of the simulated CDM halos fall with radius 
as r~ 2 . This is density profile expected for a flat rotation 
curve body. 

The competition between the dark matter and the stel- 
lar mass components within the optical radius difficults the 
analysis of solid-body portion of the rotation curve. Several 
works have argued that dark matter in dwarf galaxies domi- 
nates the total mass density profile even within their optical 
radius (Carignan & Beaulieu 1989; Broeils 1992), showing 
a constant density dark matter core (Moore 1994; Flores & 
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Figure 1. Long-slit positions superimposed on the Mrk86 B-band 
image. Close to the region B, size and position of the Mrk86-B 
bubble Ha lobes are shown (see Figures [n]a & [loja) . Relative 
coordinates (in arcsec) are referred to the r-band outer isophotes 
centre. 



Primack 1994; Salucci fc Pe rsic 1997). On the other hand, 
Lo, Sarge nt fc Young and Staveley-Smith, Davies & 

Kinman ( |l992| ) deduced reasonable virial mass to blue light 
ratios, lower than 7Mrri/ L0, fo r a large fraction of their 
samples. Loose & Thuan ( 1986 ) found that the virial mass 
of Haro 2 can be reproduced just adding the s tellar and HI 
mass components. Also the study of Swaters (1998) of the 
rotation curves of 44 dwarf galaxies indicates that the mass 
of a large fraction of these galaxies could be dominated by 
the stellar component, even at distances larger than three 
disk scale lengths. 

One of the main sources of uncertainty in all these stud- 
ies is the mass-to-light ratio adopted for the stellar com- 
ponent (Meurer, Staveley-Smith & Killeen 1998; Swarters 
1998). Therefore, high quality optical and near-infrared 
imaging and spectroscopy in order to obtain the physical 
parameters of the stellar populations and derive reasonable 
mass-to-light ratios is mandatory to prevent this inconve- 
nient which is inherent to this kind of kinematical studies. 

Superimposed on the regular solid-body portion of the 
velocity field, peculiar motions of the ionized gas have been 
observed in many star forming dwarf galaxies (Tomita et al. 
1997; Petrosian et al. 1997). They have b een commonly ex - 
plained as infalling motion s of HII regions ( Saito et al. 1992), 
multiple clouds merging ( |Skillman fc Kennicutt 1993Q and 
local p eculiar gas motions in duced by violent star formation 
events (Petrosian et al. 1997). Very high star formation rates 



associated with these intense star forming events have been 
demonstrated to be able to produce a cavity of shock-heated 
gas due to the energy input provided by supernovae and 
stellar winds (Chevalier & Clegg 1985; Vader 1986, 1987). 
This hot gas will accelerate the ambient interstellar medium 
resulting in a collective supernova-driven wind. In fact, sev- 
eral galactic supernova-driven winds phenomenae have been 
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Figure 2. Interpolated 2D velocity field. Only regions with 
Sjla >1.5x 10 — 16 ergs -1 cm -2 arcsec -2 are shown. 



found to be associated with violent star formation places in 
dwarf galaxies (Roy et al. 1991; Izotov et al. 1996; Martin 
1996, 1998, CM98 hereafter). They have been generally de- 
tected as holes in the neutral hydrogen distribution (Puche 
et al. 1992; Brinks 1994), bubbles or shells in Ha emission 
(Marlowe et al. 1995, MHW hereafter) or fro m their hot gas 
X-ray emission (Bomans, Chu & Hopp 1997). 

The existence of these phenomenae could produce 
the loss of a significant fraction of the galaxy interstellar 
medium. Depending on the final destination of the acceler- 
ated gas, these structures could produce no mass loss, blow- 
out, only affecting the galaxy chemical evolution, or blow- 
away processes, with a significant loss of interstellar mass 
(Young & Gallagher 1990; CM98; Mac Low & Ferrara 1998). 
Consequently, these supernova-driven galactic winds are ac- 
cepted to be a key parameter in the dwarf galaxy formation 
(Silk et al. 1987; Mori et al. 1997) and evolution (MHW; 
Mac Low & Ferrara 1998). 

Blue Compact Dwarf galaxies, with intense recent or on- 
going star forming activity, are those systems where the in- 
terplay between star formation and the interstellar medium 
is more feasible to be studied. However, although the ma- 
jority of the BCD galaxies are iE type BCDs (~70 per 
cent; Thuan 1991), with star formation spreading over the 
whole galaxy, the effects of the supernova- driven winds have 
been mainly studied in dwarf amorphous galaxies (see, e.g. 
MHW), which show nuclear star forming activity. 

The galaxy Mrk86=NGC 2537 (Shapl ey fc Ame s 1932; 
Markarian 1969), also known as Arp 6 (Axp 1966), con- 



stitutes an excellent laboratory to test the properties and 
effects of the supernova-driven winds on the interstellar 
medium of dwarf galaxies, as a nearby prototype of the iE 
BCD galaxies class. 

After introducing Mrk86 in Sect. ^[ we describe the 
observations and data reduction in Sect. ^. Results on 
the global velocity field of Mrk86 are given in Sect. [|. In 
Sect, [sj] we describe the evolutionary synthesis models ap- 
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Figure 3. Heliocentric velocity profiles compared with the iC-band and EWhc« profiles: a) Slit #9R. b) Slit #7R. c) Slit #2b. d) 
Slit #2r. e) Slit #4b. f ) Slit #4r. Special symbols used in panels e and f correspond to those regions marked in the left panel of the 
Figure |uj. Open circles are velocity data obtained from just one emission line, meanwhile filled circles are those measured using several 
emission lines. The peak in the EWjja profile of the panel a) has been taken as reference for relative distances in this figure and it is 
due to the contribution from the close knot #16 (G99). The same for the panels b, c and d and the knot #15 (see Figure pi). Distances 
in the panels e and f are referred to the slit region closer to the field star position (see Figure [lo[) . Lower panels show the residuals from 
the global velocity gradient fitted. 



plied. Then, we show the physical p rope rties of the Mrk86- 
A (CM98) expanding bubble (Sect. 5.2), an d th e new bub- 
bles detected Mrk86-B and Mrk86-C (Sect. EJ & ph. The 



velocity dispersion measured in Mrk86-C is also discussed 
m Sect. |. Finally, summary and conclusions are given in 
Sect. [?]. We have used ffo=50kms -1 Mpc -1 and go=0.5 
through this paper. 



2 MRK86 

The galaxy Mrk86 {a (1950)= 08 ,l 09 m 43 s ; S (1950) = 
+46° 08' 33"), as a prototype of the iE type BCD galaxies, 
is characterized by a smooth elliptical low surface bright- 
ness underlying stellar component on which several irregular 
knots of star formation are superimposed. See Figure jj] for a 
B-band image. Up to 50 different knots have been detected 
(Gil de Paz et al. 1999, G99 hereafter). In the Figure | we 
show those line-emitting regions with Ha surface brightness 
higher than 1.5xl0 -16 ergs -1 cm -2 . Assuming a receding 
velocity of 443 kms -1 (de Vaucouleurs & Pence 1980), the 



distance to Mrk86 would be 8.9 Mpc and the physical scale 
42.9 pc/". Mrk86 has an absolute magnitude of M s =-17 m , 
and it has been detected on 1.2, 2.8 and 6.3 cm (Klein, 
Wielebinski & Thuan 1984, Klein, Weiland & Brinks 1991); 
the galaxy is a source of HI 21 cm emission (Thuan & Mar- 
tin 1981; Bottinelli et al. 1984; WHISP survey, Kamphuis, 
Sijbring & van Albada 1996); also it was detected in CO 
(Verter 1985; Sage et al. 1992); it is a strong IRAS source 
(Lonsdale et al. 1985; Dultzin-Hacyan, Masegosa & Moles 



1990); near- infrared data were obtained by Thuan (1983) 
optical images were obtained by Hodge fc Kennicutt (1983) 



Loose & Thuan (1985) and Ojha & Joshi (1991); and, finally, 



ultraviolet spect ra has been analy zed by Fanelli, O'Connell 
& Thuan (|l988[) and Longo et al. (|l99l|). 



Papaderos et al. (1996a, 1996b; P96a & P96b hereafter) 
obtained the physical parameters of the plateau, central 
and exponential photometric components. These compo- 
nents have been related with three different stellar popu- 
lations (respectively, currently forming, ~lGyr and >7Gyr 
old; Gil de Paz, Zamorano & Gallego 1998, G98 hereafter). 
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Table 1. Journal of observations. 







Spectroscopic observations 






Exp. 


Slit 


Spectral 


Disper. 


Telescope 


time (s) 




range(nm) 


(A/pix.) 


CAHA 2.2-m 


3600 


1,2,4,6b 


330-580 


2.65 


CAHA 2.2-m 


1800 


5b 


330-580 


2.65 


CAHA 2.2-m 


3600 


l,2,4,5,6r 


435-704 


2.65 


CAHA 2.2-m 


3600 


3 


390-650 


2.65 


INT 2.5-m 


1800 


7,8R 


637-677 


0.39 


INT 2.5-m 


900 


9R 


637-677 


0.39 






Image observations 






Exp 


Filter 


Scale 


PSF 


Telescope 


time (s) 




("/pixel) 


(arcsec) 


JKT 1.0-m 


600 


B 


CK'33 


l'.'O 


CAHA 1.5-m 


2400 


V 


rt'33 


1'.'6 


CAHA 2.2-m 


600 


Gunn-r 


C/27 


l'/9 


INT 2.5-m 


900 


Ha 


OC'57 


2'.'3 


INT 2.5-m 


900 


[OIII] 


0757 


2'.' 5 


UKIRT 3.8-m 


100 


K 


0'.'60 


1'.'4 



Table 2. Velocity gradients. u;o are the velocity gradients mea- 
sured along the slits, meanwhile u> are the values projected to 
PA=12°. Distances along the rotation axis, d, have been mea- 
sured relative to the r-band outer isophotes centre. The central 
value has been estimated from the interpolated 2D velocity map. 



3 OBSERVATIONS AND REDUCTIONS 

We have obtained a total of 14 long-slit optical spectra in 

Jour- 
Low- 



10 different orientations and positions (see Figure h j 
nal of spectroscopic observations is given in Table 1 
intermediate resolution spectra (FWHM=6A in the light of 
Ha) with dispersion 2.65 A/pixel were obtained with the 
Boiler & Chivens spectrograph at the Cassegrain focus of 
the 2.2-m telescope at Calar Alto (Almerfa, Spain) in Jan- 
uary 1993. We used a 1024x1024 Tek 24^m CCD. High 
resolution spectra (FWHM=0.9A in Ha and 0.39 A/pixel; 
slits #7R, #8R and #9R) were obtained using the IDS in- 
strument at the Isaac Newton Telescope (INT) of the Ob- 
servatorio del Roque de los Muchachos (La Palma, Spain) 
in January 1998. The detector was a 1024x1024 Tek 24/xm 
CCD. The spectra were reduced, wavelength and flux cali- 
brated, using standard FIGARO (January 1993) and IRAF^ 
procedures (January 1998), and making use of the corre- 
sponding standard stars and lamp calibration spectra. Line 
fluxes were measured using the IRAF task SPLOT, obtaining 
errors lower than 15 per cent. 

We have also obtained B, V, r and A'-band images in 
different observing runs as described in Table |l|. Gunn-r 
band image was obtained at the 2.2-m telescope at Calar 
Alto in February 1992 with a GEC CCD with 800x1156 
30/j.m pixels. The Johnson-!/ band image was obtained at 
the 1.52m telescope at Calar Alto in December 1993 with a 
1024x1024 19/xm pixels CCD. B band image was obtained 
at the 1-m Jacobus Kapteyn Telescope (JKT) at the Ob- 
servatorio del Roque de los Muchachos in November 1997 
with a 24^m 1024x1024 pixels Tek CCD. Finally, the near- 
infrared if-band image was obtained at the 3.8-m UKIRT 
telescope at Mauna Kea observatory in April 1993. We used 
the 62x58 InSb array camera IRC AM. Broad band im- 



* IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which are operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 



d (kpc) 


t^o (km s 1 kpc 1 ) 
Steep/flat comp. 


u) (km s 1 kpc 1 ) 
Steep/flat comp. 


Slit 


1.03 


55±4/8±3 


56±5/8±3 


#9R 


0.85 


69±5/18±2 


70±5/19±2 


#7R 


0.73 


21±4 


21±4 


#2r 


0.73 


22±5 


22±5 


#2b 


0.00 




34±5 




-0.12 


33±4 


33±4 


#4r 



ages were reduced using standard MIDAS and IRAF pro- 
cedures. The procedures described by Aragon-Salamanca et 
al. (1993) were employed for the reduction of the Jf-band 



near-infrared image. All the broad-band images were cali- 
brated observing standard stars at different airmasses. 

Additional narrow-band images were secured for us dur- 
ing service time at the 2.5-m Isaac Newton Telescope at the 
Observatorio del Roque de los Muchachos in December 1993. 
The detector was an EEV CCD with 1280x1180 22.5^m 
pixels. [OIII]5007A (A =5014A; FWHM=50A) and on/off 
Ha (Ag n =6556A, FWHM=60A; Ajf =6607A, FWHM=53A) 
images were obtained. The continuum subtraction was per- 
formed using, respectively for [OIII] and Ha, the V-band 
and off-Ha images. 

The [OIII] and Ha images were flux calibrated as fol- 
lows. First, we convolved the b {blue-arm; see Table [j]) 
spectra with the [OIII] filter and the r {red-arm) spectra 
with the Ha filter (see Table [j] & Figure |). Then, we col- 
lapsed these spectra in wavelength, obtaining the flux cal- 
ibrated spatial profile of the [OIII] and Ha emission. Now, 
we cropped and averaged in the [OIII] and Ha images those 
regions covered by the b and r slits, obtaining the [OIII] and 
Ha spatial profiles given by the images. Comparing these 
spatial profiles with those obtained from the spectra we cal- 
ibrated both the [OIII] and Ha images. Finally, these cal- 
ibration relations were corrected for the sensitivity of the 
filters at the corresponding wavelength and, in the case of 
the Ha emission, from the contribution of the [NII]6548A 
and 6583A lines. The reliability of this method was demon- 
strated after applying it to all the blue-arm and red-arm 
spectra, obtaining similar results within an error of 10 per 
cent. 



4 GLOBAL VELOCITY FIELD 

We have obtained the spatial variation of the ionized 
gas velocities using the long-slit spectra. We have used 
the rvidlines IRAF task measuring the [OII]3727A, H/3, 
[OIIIJ4959A and [OIII]5007A line velocities for the blue- 
arm spectra, and [OIII]4959A, [OIII]5007A, [NII]6548A, 
Ha, [NII]6583A, [SII]6717A and [SII]673lA for the red-arm 
spectra, weighted with their relative intensities. No signifi- 
cant differences were observed using allowed and forbidden 
lines separately. Errors in the velocity were estimated there 
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Figure 4. Thin lines: Neutral hydrogen (dashed), stellar- 
underlying (dotted), stellar-starburst (dot-dot-dot-dashed), dark 
matter (dot-dashed) and total (solid) mass density profiles for 
i?DM=0.5 kpc. Thick lines: Dark matter (dot-dashed) and total 
(solid) mass density profiles for Rr)M=3.0kpc. 
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Figure 5. Position angles, i and i c , and coordinates x and y for 
two equatorial HII regions (#1 & #2). Solid plane represents the 
galactic equator. 



mass density distribution, related with the presence of these 
massive star forming regions. However, a merging with an- 
other dwarf galaxy or gas cloud with an independent velocity 
field, as it has been proposed to explai n the velocity fie ld of 
I Zw 18 dSkillman fc Kennicutt 1993|) and II Zw 40 ( [Van 
Zee et al. 1995), should not be ruled out. In that case, this 



where several emission lines could be measured (see Fig- 
ure |). 

From the velocities determined in the spectra, and the 
slit positions given in Figurejll, we have reconstructed a 
2D velocity map (see Figure fusing the IRAF task xyz- 
toim. A careful examination of the interpolated 2D ve- 
locity field yield a global solid-body velocity field of maxi- 
mum velocity gradient f2=34kms _1 kpc -1 and orientation 
PA~12° (rotation axis PA~ -78°). The heliocentric ve- 
locity measured at the galactic centre, as given by the r- 
band outer isophotes centre , is a bout 470 km s -1 . Although 
De Vaoucoleurs & Pence (1980) gave a median heliocen- 
tric velocity of 443kms _1 , a detailed examination of their 
Fabry-Perot interferogram shows that the velocity close 
to the galactic center is significantly higher, between 450- 
460 kms" 1 . 

In Figure ^ we show the velocity profiles measured along 
the slits #9R, #7R, #2b, #2r, #4b and #4r, those posi- 
tioned closer to PA=12°. The six panels have been arranged 
in slit position order, from West to East, where slit #4b is 
that closely crosses the galactic centre, as defined by the r- 
band outer isophotes (a (1950)= 08 h 09 m 42=56; 8 (1950) = 
+46° 08' 33'.' 8). Northernmost part of the slit is represented 
at the leftside of each panel. 

The panels a) and b) (slits #9R and #7R) show 
two different velocity components, one very steep gradient 
and a second flatter velocity gradient. The steeper com- 
ponents seem to be associated with regions of enhanced 
if-band luminosity and relatively high star forming activ- 
ity (EWhq ~100A), the emission knots #15 [a (1950) = 
08 h 09 m 41=17; S (1950)= +46° 08' 45'.' 1; G99) and #16 (a 
(1950)= 08 h 09 m 41 s 00; S (1950)= +46° 08' 43'.'0; G99; see 
Figure ^) . This increase in the radial velocity gradient could 
be due to an enhancement of the mass density over the global 



merging could be responsible for the triggering of the star 
formation in the emission knots #15 and #16 (G99). 

On the other hand, the #9R and #7R flat components 
and the gradients measured along the slits #2b, #2r and 
#4r should be related with a more relaxed global veloc- 
ity field. In Table ^ we give these velocity gradients, being 
ujo the gradients measured along the slits and ui the values 
projected to PA=12°, i.e. uio=^> x cos(PA— 12°). Velocities 
affected by local motions, such as supernova-driven winds 
(see Sect. p|) have not been included in the calculus of the 
velocity gradients. From Table ^ we see that the flat compo- 
nent of the velocity gradient grows from 10-20 km s" 1 kpc" 1 
in the outer galaxy regions to ~34kms _1 kpc -1 close to the 
galactic centre. 



4.1 Mass density profile 

We will assume that the motion of the HII regions which 
produce the line-emission observed is due to rotation. The 
slit #4r was placed very close (~3") along the galactic equa- 
tor. Therefore, we could compare its velocity profil e with the 
radial compone nt of the circular velocity curve (Binney & 
Tremaine 1987). We will parameterize the mass distribu- 



tion of the galaxy (stellar and dark matter, and neutral and 
molecular hydrogen). We have assumed that the photomet- 
ric centre, given by the r-band outer isophotes, coincides 
with the kinematical centre (see Ostlin et al. 1998). In order 
to compare the velocity curve obtained along the slit #4r 
with the radial component of the circular velocity curve we 
will adopt the radial velocity measured close to the galactic 
center as the sistemic velocity. 

Stellar. From the analysis of the surface brightness pro- 
files of 14 Blue Compact Dwarf galaxies, P96a found that 
these profiles can be fitted using three distinct components, 
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Figure 6. Total (solid) and stellar (dotted) radial velocity curve 
predicted for i c =0° and inclination of the rotations axis, i, 0, 
30 and 60°. Black and white dots are, respectively, receding and 
approaching velocities. The dot-dot-dot-dashed, dot-dashed and 
dashed lines, respectively the starburst, DM and HI density com- 
ponents are represented just for i=0° and i?DM=ikpc. 



the underlying, plateau and starburst components. Using a 
numerical deprojection procedure, they obtained the cen- 
tral luminosity densities and scale lengths for these com- 
ponents. G98 estimated the ages for the starburst and un- 
derlying components to be, respectively, 1 and 7Gyr, with 
one fifth solar metallicity. Then, we have taken the mass- 



to- light ratios predicted by the Bruzual & Chariot ( 199C ) 
evolutionary synthesis models for these ages and metallici- 
ties, under the assumptio n of instant aneous star formation 
and using an Scalo IMF (3calo 1986). These mass-to-light 
ratios result in M/Ls=0.5 Mq/Lq for the starburst compo- 
nent and M/L_b=2.6 Mq/Lq for the underlying component. 
Although quite significant in luminosity, the plateau compo- 
nent represents a very small fraction in mass because its lu- 
minosity is produced by regions with star formation more re- 
cent than 20 Myr (G9 9) and mass-to-light ratios lower than 
M/L s =0. 08 Mq/Lq flBruzual fc Chariot 199^ ). 

Finally, using the parameters given by P96a for the 
Mrk86 starburst and underlying components we obtain 



under / \ 

p M 



under , 

Po exp ( 



sb / \ 

p i r ) 



sb , 

Po exp( 



-r/Rs 



(1) 

(2) 



-3 



where 

under on ~± ~™ 

Po —20 atn cm 

Runder = 1.0 kpC 

po 5 =31 at H cm -3 
7? s6 =0.14kpc 



Neutral hydrogen. If we assume that the HI distribution, 
as given by the WHISP survey^, falls approximately as an 



t 



http: //thales . astro . rug.nl/~whisp/Database/ 



Figure 7. Radial velocity curves predicted for i=0° and 
ic=0,30,60°. The meaning of the symbols and line patterns are 
the same that in Figure pi Starburst, DM and HI density compo- 
nents are represented just for i e =0° and -Ro.M = lkpc. 



exponential function, and using the same deprojection pro- 
cedure that for the optical data (see P96a), we obtain 



HI/ \ HI / /r> \ 

P (r) = Po exp(-r/Rm) 



(3) 



where 
,R H i^2.0kpc 



pJp~0.7at H cm -3 



Molecular hydrogen. Sage et al. (1992) gave a 
M(H 2 )/M(HI) ratio for Mrk86 of 0.03. Therefore, we can 
assume negligible the effect of the molecular hydrogen in 
the gravitational potential of Mrk86. 

Dark matter. Finally, we have considered the effect of the 
dark matter distribution over the global velocity field. Fol- 
lowing the most recent observational works (Salucci & Persic 
1997; Flores & Primack 1994; Moore 1994) and the study 



of Navarro et al. (1996) concerning the effects of supernova- 



driven winds over standard CDM profiles, we have assumed 
the existence of a central core in our dark matter density 
profile. Both the density distribution of a modified isother- 
mal sphere (see Binney fc Trem aine 1987) and the universal 
profile given by Burkert ( 1995 ) include this central core and 
reproduce reasonably the observations. We have used for this 
work the simplier mo dified isothermal profile, also following 
Mac Low & Ferrara Jl998|), 



DM 

DM, % Po 
P M = 



1 + (r/RoM) 2 



(4) 



Burkert (1995; see also Mac Low & Ferrara 1998) has shown 
that the central density pp is related with the scale radius 
Rdm through the expression 



DM 

Po 



2.7 x 10 



Rt 



kpc 



-2/3 



M kpc 



(5) 
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Figure 8. Change in the radial velocity gradients with the dis- 
tance a 
Table 



ong the rotation axis. Data points have been taken from 
. For d=0 we have taken the velocity gradient measured 
from the interpolated 2D velocity map. Thick-lines represent the 
velocity gradient predictions for 30° (solid), 45° (dashed) and 
60° (dot-das hed) inclined rotation axis with the density profile 
given in Sect. [f. l| . Thin-lines are the predictions for 2=30° (solid) 
and 45° (dashed), but using i? unt j er =0.5 kpc, instead of 1.0 kpc. 
White dot was obtained for a projected distance of —0.12 kpc east 
of the galactic centre (see Figure [j]). 



Thus, this profile has only one free parameter, i?DM- 

In Figure ^ we show the total and DM mass density 
profiles for J?DM=0.5kpc and 3 kpc. In both cases we see 
that the total mass density profile is dominated by the stellar 
component within the central 4 kpc. 

Now, we could compare the radial component of the 
circular velocity curve derived from the total density profile 
with the radial velocity data obtained along the slit #4r, 
placed close along the galactic equator. However, since the 
ionized gas emission is produced in several individual HII 
regions, the projected radial component of the modelized 
circular velocity will depend on the position of the corre- 
sponding HII region. 



4.2 Ionized gas geometry 

We will compare the velocity data obtained along the equa- 
torial slit #4r with the projection of the rotation curve 
adopting several geometries. 

We will assume that the ionized gas is distributed in a 
thin-disk with inclination i. This angle measures the inclina- 
tion of its rotation axis with regard to the plane of the sky 
(i.e. i—0° for a edge-on disk; see Figure ^). In the Figure || 
we also show how the angles i and i e and the coordinates 
x and y describe the position in the disk of an HII region. 
Under the hypothesis of thin-disk, if i 7^0° only those regions 
with angle i e ~0° will be observed through the slit #4r. This 
is the case shown in the Figure |(| 

However, if the galaxy would have an inclination close to 
i=0°, HII regions with different position angle, i e , could be 
observed through the slit. In this situation, we could observe 
regions with different radial velocity at the same apparent 



position (regions #1 & #2 in Figure |5|). For example, if we 
consider the mass density profile given in Sect. [ll] (with 
flj3M = lkpc), we could measure changes in the radial ve- 
locity of Wkms" 1 between two regions with x=0kpc and 
a;=3kpc for i/=2.5kpc (see Figure m. Therefore, the exis- 
tence of inhomogeneities in the ionized gas distribution in 
the galaxy could also produce short scale variations in its 
velocity field. 

The emitting HII regions can have, in principle, any 
position angle ji c | <180°. In fact, a region with a given i e at 
present, would evolve toward larger i c angles in the future. 
Had all the HII regions be placed in planes of constant angle 
i e , the resulting rotation curves would be the ones shown in 
Figure f?| 

From the Figure ^ we see that the stellar mass density 
component dominates at the inner 4 kpc. However, the radial 
velocity data can be fitted adopting different geometries for 
the ionized gas distribution (see Figures |and0). We could 
reproduce the velocities measured if z=0° and 30° < i c <70°, 
or if i~40°and i c =0° . 

On the other hand, the way the velocity gradient per- 
pendicular to the rotation axis, u>, decrease from the equator 
to galactic outer regions (see Table |^) resembles the velocity 
fields observed in spiral galaxies with intermediate inclina- 
tion (e.g. Giovanelli & Haynes 1988). We will test this point 
assuming that the emission observed effectively comes from 
a rotating thin-disk with gravitational potential tha t given 
by the mass distribution described before (Sect. 4.1). Then, 



using this gravitational potential we have estimated the ex- 
pected radial velocity gradients for different distances, d, 
along the rotation axis (see Table [2J). In Figure |§| we show 
the gradients expected, measured as the change in radial ve- 
locity of their inner two kiloparsecs. Data points have been 
taken from Table ^| The value obtained for the slit #9R is 
very uncertain since the small number of points used to fit 
its velocity gradient (see Figure ^). The zero value repre- 
sents the gradient along the galactic equator obtained from 
the interpolated 2D velocity map. From this figure we note 
that, although a decreasing in the velocity gradient is effec- 
tively observed, we cannot reproduce the values measured. 
The more feasible explanation is that the ionized gas is prob- 
ably distributed in a relatively thick-disk. However, the use 
of a steeper mass density profile could also reproduce the 
velocity gradients measured (see Figure ^). 

Therefore, we conclude that the observed global velocity 
field of Mrk86 can be reproduced if the stellar component 
dominates the total mass profile within its optical radius, 
and if the emitting ionized gas is distributed in a probably 
thick inclined disk. The inclination of the disk with regard 
to the plane of the sky will be about 50° (i~40°). 



5 KINEMATIC AL EVIDENCES FOR 
SUPERNOVA-DRIVEN WINDS 

Several under-kiloparsec-scale kinematical structures are ob- 
served superimposed on the global velocity gradients pre- 
viously described. These kind of kinematical disturbances 
have been widely observed in dwarf galaxies (see Tomita et 
al. 1997; Petrosian et al. 1997). Different mechanisms could 
properly explain thes e disturbances, in cluding infalling mo- 
tions of HII regions ([Saito et al. 1992 ) or inhomogeneities 
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in the ionized gas distribution (see Sect. kf). However, in our 
case, the deep minima (see Figures ^b, ge, ^f, [l2^ and |l2]b) 
and velocity maximum (Figure ^) observed are clearly spa- 
tially correlated with intense star forming knots, in dica t ing 
th e ex istence of supernova-driven winds (see Sect. 



5.4) 



5.2 



5.3 



5.1 Description of the evolutionary synthesis 
models 

In order to determine the physical properties of the 
star forming knots, in particular, those associated with 
supernova-driven wind phenomenae, we have made use of 
the Leitherer and Heckman (1995; LH hereafter) models. 
These single burst models do not take into account any con- 
tribution from the old underlying population neither that 
from emission lines (i.e. [OIII]5007A, Ha). 

Optical and near infrared colours measured in the 
outer galaxy regions (B - V=0.8, V - r=0.0, V - K=2.5; 
G98), where the exponential underlying population domi- 
nates (P96a), allow to estimate an age of about 7Gyr (G98 
using the Bruzual & Chariot 1996 models) for this under- 
lying component. Assuming different burst-strength values, 
ranging from 0.001 (0.1 per cent in mass) to 1, we have re- 
built a complete set of LH instantaneous burst models o f 
0.25 and 0.1 Zq metallicities, Salpeter IMF ( fcalpeter 195. 



with m= -2.35, M; o iu = 1Mq and M up = 100M Q . Then, 



relation ( |Kent 1985| ), 



r = R + 0.41 + 0.21 x(V-R) 



the 



(6) 



has been applied in order to transform from Johnson- R to 
Gunn-r magnitudes. Finally, using the Ha equivalent widths 
measured from the spectra, we have corrected the Gunn-r 
magnitudes from Ha emission. 



5.2 Previously detected Mrk86~A bubble 

Using long-slit echelle spectra in Ha+[NII] with llkms -1 
resolution (FWHM), CM98 inferred the presence of an 
expanding bubble extending 22" (940 pc) westward of the 
galaxy west end. The physical parameters of this bubble are 
given in Table [| In Figure ^| we show the A-big bubble re- 
gion in Ha and S-band. Following CM98, we have taken the 
knot marked in Figure ^ as A s b (#31 emission knot in G99) 
as the starburst precursor. 

For this region and for the Mrk86-B and C bubble 
precursors, we have determined physical apertures using 
the program COBRA (G99). This program determines and 
subtracts the underlying emission. The apertures include 
those pixels with emission more intense than 1/e times 
the emission maximum. Optical-nIR colours have been ob- 
tained for these apertures and corrected for inter nal ex - 
tinction applying the relation given by Calzetti (1997), 
E(B - V) C ontinuum=0.44 xE(B - V) gas . Using the Ha-H/3 
Balmer decrements and assuming a Mathis (1990) extinction 
law, colour excesses E(B — V) gas have been obtained. The 



corrected colours are given in Table |§| These colours allow 
to estimate the mean physical properties, i.e. age, metallic- 
ity and burst strength for these regions. In Table j3| we show 
the physical parameters of the LH model that better fits the 
colours measured. 




10 -10 -20 
Relative Distance (") 

Figure 9. Mrk86-A region. B image and Ha contours are shown. 
Displacements between S-band knots and Ha maxima are evi- 
dent in the position of the Mrk86-A bubble starburst precursor 
(CM98), A sb . The limits of the Mrk86-A Ha emission given by 
CM98 are also shown, joint with the positions of the #15 and 
#16 (G99) emission knots. 



In the Mrk86-A case, the colours measured yield a LH 
model with 0.25 Zq metallicity, burst strenght of 1 per cent 
and age of about 14 Myr (see Table |^) . This value is slightly 
different from the dynamical time deduced by CM98 (con- 
verted to #o=50 kms -1 Mpc" 1 ), that is t dyn ^12Myr. This 
discrepancy could be explained taking into account that 
tdyn describes the age of the expanding bubble, that dif- 
fers from the starburst age in a time similar to the more 
massive stars main sequence phase duration. After the red 
super giant phase has started (~4x 10 6 yr after the starburst 
for M~40 Mq; Maeder 1990), massive stars winds and later 
supernova explosions would commence to take place and a 
cavity of shock-heated gas could begin to form. 



5.3 Mrk86-B 

A local velocity maximum and minimum is observed along 
the #4b slit velocity profile (see Figure |j|e). They differ 
68kms _1 in velocity and are 10" (430 pc) apart. These fea- 
tures are compatible with the approaching (north lobe) and 
receding parts (south lobe) of a bipolar expanding struc- 
ture with a characteristic projected expansion velocity of 
34kms _1 . The approaching velocity minimum is also ob- 
served in the #4r slit velocity profile (see Figure |^f). As- 
sociated with these kinematical features a clear bubble-like 
structure is observed in the Ha narrow band image (see Fig- 
ure [l(]). Assuming a distance of 8.9 Mpc, a physical size of 
750x510 pc 2 is obtained. This structure resembles those ob- 



served in I Zw 18 (Martin 1996), M82 (Heckman et al. 199C) 



and NGC 1705 (MHW). The difference in size between both 
lobes could be related with a strong ambient density gradi- 
ent, with higher density toward the south lobe region. 

At 3" (130 pc) to the NE from the bubble geometrical 
centre (origin of coordinates in Figure |lo| ) an optical emis- 
sion knot (#49 in G99) is observed in the B, V, r and K 
bands (Figure [h], left panel). If we compare the instanta- 
neous burst evolutionary synthesis LH models to the colours 
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Figure 10. Left panel: Mrk86-B region in B-band (image) and Ha (contours). Region 1" East and 9" North is a field star. Right panel: 
Ha (image) and [OIII]/Ha (contours). See Figures and |jf for the velocities measured in the regions marked on the left panel. In the 
[OIII]/Ha case, some contour-map highly noisy regions have been artificially removed. 



Table 3. Bubbles physical parameters. dE/dtki n values given in the upper panel are deduced applying the expressions 
of rabble and i>exp from Castor, McCray & Weaver ( |1975[ ). Lower panel shows the parameters of the best fitted 
evolutionary LH model and the predicted collisionally excited Ha luminosity. Optical-nIR colours in parenthesis 
( upper vanel) a re those given by the corresponding LH model. E(S — V) C ontinuum is assumed to be 0.44xE(B — V) gSjB 
(Oalzctti 1997). E(B — V) gas value s for Mrk86— A bubble have been obtained by CM98 and this work (in parenthesis). 
Following the Castor et al. (1975) bubble model, dynamical ages, tdym have been estimated as 0.6xrb u bbleAfexp 
(see, e.g. CM98). 
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Figure 11. Mrk86-B region, a) Filled dots: #4b long-slit velocity profile. Upper panel: Observed [OII]/[OIII] profile (solid line), [OIIIJ/H/3 
profile (dashed line). Lower panel: [Oil] (solid), H/3 (dashed) and [OIII] (dot-dashed) profiles, b) Filled dots: #4r long-slit velocity profile. 
Upper panel: Observed [OIII]/Ha profile (solid), [SII]/Ho profile (dashed), [Nil] /Ho profile (dot-dashed) [SII]6717A/[SII]673lA profile 
(dotted). Lower panel: H/3 (solid), [OIII] (dashed), Ha (dot-dashed), [Nil] (dotted), [SII] (thick solid). Relative distances are referred to 
the same regions as in Figures He & Hf. 



measured (see Sect. 5.1 and Table |3|), considering an age for 
the old underlying population of 7Gyr (G98), the best fit 
obtained yields an starburst age of about 11.5 Myr, a burst 
strength of 1 per cent and a metallicity of about 0.1 Zq. We 
are confident with the instantaneous burst assumption be- 
cause we are dealing with single star forming regions (LH). 
From the B-band total luminosity, and taking into account 
the burst strength derived, we estimate the involved total 
burst mass to be about 6.3 xlO 4 Mq. 

Assuming that the kinematical features observed rep- 
resent the foreground and receding parts of the expand- 
ing bubble, the dynamical age for this bubble, obtained as 
£dyn=0.6x(rs + ri<t)/(vs +Un)j will be about 6.5 Myr, quite 
similar to the evolution time deduced for the starburst re- 
gion, if we consider the massive stars main sequence time, 
i.e. about 4-5 Myr, subtracting it from the starburst age 
(11.5 Myr). 

Generally, the Ha luminosity observed in expanding 
bubbles is accepted to be produced by photoionization from 
massive stars of the starburst precursor (MHW; Lehnert & 
Heckman 1996; Martin 1997). Then, the predicted collision- 
ally excited Ha contribution should be irrelevant to the total 
Ha emission in the expanding lobes. The Ha luminosities 
measured for both lobes are LH° rth =4.7x 10 37 erg s -1 and 
L|^ th =2.5xl0 38 ergs -1 . Considering that both lobes are 
well reproduced by ellipsoids of revolution, we can obtain 
the lobe surface areas, that allow to estimate the colission- 
ally produced Ha luminosity 



where no is the HI ambient density, t» exp the bubble ex- 
pansion velocity, / is the number of Ha photons produced 
per shocked proton, and un a is the frequency of the Ha 
photons. The predicted collisionally excited luminosities are 



1 =2.4xlO JD ergs~ 



and L|° a uth =1.6xl0 3 



ergs 



less 



Lhc — Slobe n V exp f H l>Ha 



(7) 



than the 10 and 1 per cent of the total Ha luminosity. We 
have assumed that approximately 0.1 Ha photons are pro- 
duced per shocked proto n. This is the value predicted by 
the Shull & McKee ( J197S| ) models for a 50kms _1 shock. We 
have also adopted a HI ambient density of no ~0.3 cm -3 , 
that is the same value adopted by MHW for a sample of 
dwarf amorphous galaxies. This value is obtained consid- 
ering the galactocentric distance of the Mrk86-B bubble, 
1.2 kpc, and the galaxy HI density profile given by the equa- 
tion pi 

Due to the uncertainties present in the HI total mass 
and scale measurements, we could be dealing with signifi- 
cantly higher HI ambient densities (e.g. no ~lcm -3 ), that 
could enhance the predicted collisional Ha luminosity. In 
any case, this contribution would not dominate the total 
Ha luminosity in the north lobe (< 30%), and would be 
negligible in the south lobe case (~ 3%). 

Then, if the emission is photoionization dominated, the 
ratio between the luminosity of both lobes has to be re- 
lated with the solid angle subtended from the starburst by 
each lobe. In fact, assuming revolution symmetry for both 
lobes, the solid angle ratio obtained is fisouth/fiNorth=5.5 
(^South=6.6 sr and fiNorth — 1-2 sr), very close to the ratio 
■ South /T Norths 3 In the rigM panel of Figure [Jol we show 



IT No 
/ L Ha 
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Figure 12. Mrk86-C region, a) ^£7R and b) ^£8R long-slit velocity profiles. Uppev panel: Velocity dispersion corrected from (^instrumental • 
Lower panel: Filled circles represent the velocity data subtracted from large scale velocity gradients. Dotted circles correspond to those 
velocities measured using just one emission line. ±1<t errors bars have been plotted. Thin solid line represents the fit to an isotropic 
expanding bubble with w e xp=17.3 km s — 1 and radii 90 and 130 pc, respectively. Thick lines show the [Nil] /Ha (a) and [SII]/Ha (b) line 
ratios profiles relative to the bubble geometrical centre. 



the aperture angles assumed, and the position of the ionizing 
radiation source (#49 knot in G99). 

In order to estimate the relevance of the shock exci- 
tation mechanism, we have also studied the 
[OII]3727A/[OIII]5007A, [OIII]5007A/H/3, [NII]6583A/Ha 
and ([SII]6717A+673lA)/Ha line ratios across both lobes 
(see Figures & 0b). For the south lobe region we ob- 
tain low excitation, log([OIII]/H/3)~0.3, but with moderate 
log([OII]/[OIII])~0.3 values. These values can be explained 
using photoionization models with low ionization parame- 
ter (log U~— 4; Martin 1997). However, for the more remote 
north lobe regions the very high log([OII]/[OIII])~0.6 (Fig- 
ure |ll|a.) and high log([NII]/Ha)~-0.5 and log([SII]/Ha)~0 
(Figure [ll]b) values obtained imply that an additional exci- 
tation mechanism should be present (Lehnert & Heckman 
1996; Martin 1997). These line ratios are well reproduced 
taking into account the emission predicted from the models 
of Shull & McKee ( fl979| ) with shock velocities of 90 kms~\ 
This result is consistent with the higher collisional Ha lu- 
minosity predicted for the north lobe of the bubble. 



5.4 Mrk86-C 

From the high-resolution spectra #7R and #8R we infer the 
presence of deep velocity minima close to the G99 #9 knot 
emission maximum (C region in Figures ^ and |^) . These min- 
ima observed in both spectra #7R and #8R (see Figures [l^a 
& |l2p) are well understood if we are dealing with the fore- 



ground part of a highly extinct expanding bubble, with 
median parameters, rbubbie=H0 pc and v C xp=17.3kms _1 
(tdyn —4 Myr) . These parameters (see Table fi) are obtained 
fitting the radial velocity profiles (global velocity gradient 
subtracted; v x ) to an isotropic expanding structure velocity 
law (see Figures n3a & |l2|b) 



r 2 

' bubble 



(8) 



The central starburst should have a high extinction, large 
enough to obscure the Ha emission of the receding part of 
the bubble. In order to test this hypothesis we should es- 
timate the physical properties, extinction included, of the 
starburst which has originated this collective wind. 

The Mrk86-C bubble has similar appearance in all the 
optical and nIR images, including the Ha image (see Fig- 
ure ^ , in contrast with the B bubble case, probably due to 
its lower dynamical time. Its Ha emission is slightly more 
extended along the W-E than that observed in broad-band 
filters. The corresponding optical and nIR colours measured 
for the starburst are given in the Table ^ and are fully com- 
patible with a 7.7Myr old and Z=0.1Zq metallic region, 
adopting, in order to better fit the colours measured, an 
colour excess of E(£> — V) sas =0-73. 

From these parameters, and the total B absolute mag- 
nitude measured, applying the LH results we obtain a 
starburst mass of about 4.7x10 4 Mq. In addition, be- 
cause the total Ha luminosity measured for Mrk86-C, 
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1.7xl0 39 ergs -1 , is 6xl0 3 times higher than that predicted 
for the collisional case, 2.7xl0 35 ergs -1 , we can conclude 
that the line emission is also in Mrk86-C mainly due to 
photoionized gas. 

Finally, we have studied the [Nil] /Ha and [SII]/Ha 
line ratios measured from the high resolution spectra. They 
change between ~0.10 in regions close to the starburst 
centre, and ~0.25 and ~0.30, respectively, in the external 
low Ha surface brightness bubble regions (see Figures |l2fi 
& p^b). The increase in line ratios values towards bub- 
ble outer regions seems again to indicate a decrease in the 
ionization parameter due to the dilution of the radiation 
from the centralized source (Martin 1997). The values ob- 
tained are compatible with a pure low-metallicity HII region 
photoionization mechanism (see, e.g. Shields & Filippenko 
1990). 



6 VELOCITY DISPERSION IN MRK86-C 



about 10 J Mq. We can conclude that the burst total mass 
is significantly lower that the virialized mass deduced. The 
similarity of the velocity dispersion measured pixel by pixel 
rules out the global velocity gradient as origin for the ob- 
served velocity dispersion. First reason for Mrk86-C not to 
be bound arises from the existence of local gas motions in- 
duced by the Mrk86-C supernova-driven wind. Another ar- 
gument is Mrk86-C is not placed at the galactic centre and, 
therefore, is subject of strong tidal forces that prevents it to 
be virialized. 

We have obtained the turbulent velocity dispersion [[] 
(see, e.g. Fuentes-Masip 1997), applying the expression 



2 _ 2 

^turbulent ^observed 



■ °" the 



^instrumental 

(10) 



We have measured the ionized gas velocity dispersion tor 
the #7R and #8R high-resolution spectra. From the cal- 
ibration lamp and night sky lines we have estimated the 
spectral resolution in o-i nstrumcnta i = 16.4±2.1 kms -1 . Then, 
we have obtained the velocity dispersion of the Mrk86- 
C bubble from these spectra, subtracting to o- 2 bscrvcd the 
value of 0f nB t rum ental- The spatial variation of these veloc- 
ity dispersion measured is shown in Figure [l2| The ve- 
locity dispersion obtained for the whole (4'.'6 x 1") bubble 
C were o#7R=16±2kms -1 and a#8R=15±2kms -1 . In ad- 
dition, we have obtained the mean velocity dispersion for 
the whole emitting region, averaging the values measured 
in intervals of 0'.'33, 1" and 2" along the spatial direction. 
We have estimated 17±3 and 14.9±1.6 kms -1 , respectively 
for the #7R and #8R spectra and 0'.'33 intervals, 18±5 
and 15.4±0.5 kms -1 for 1" intervals, and finally, 19±5 and 
15.2i0.3kms- 1 for 2" intervals. 

We observe that there is no significant differences be- 
tween local and total velocity dispersion values, i.e. the ve- 
locity widths obtained for regions of size 0'.'33x 1" are similar 
to that obtained for an important fraction of the whole emit- 
ting region (4'.'6 x 1"), and that they do not change with the 
considered spatial interval. 

We have now considered the expresion for the mass of 
a virialized system given by Guzman et al. (1996, 1997; see 
also Bender, Burstein & Faber 1992 and Gallego et al. 1998), 



that yields a (Tturbulent — 10.7kms -1 and 10.0 kms -1 , re- 
spectively for #7R a nd #8R spectra , assuming for Ha an 
o"intrinsic=7.1 kms -1 ( Hippelcin 1986), and a temperature, 
T=10 4 K. Adopting a sound velocity for the ionized gas 
in H II regions of about 12.8 kms -1 (Dyson & Williams 
198C), the corresponding cr sounc i is approximately 8 kms -1 



M vir = 1.184 10 6 i? c (kpc) cr 2 (km/s) 2 M G 



© 



(9) 



If we take for Mrk86-C the S-band effective radius of 
the starburst measured with COBRA, R c ~36pc, and 
cr~15 kms -1 , we obtain a virialized mass (Eq. ^j) for Mrk86- 
C of about M~ 1O 7 M . 

Let us compare this result with the Mrk86-C starburst 
total mass estimated from the mass in newly formed stars 
(M^°" s " g =5.3x 10 4 Mq), underlying stellar population, gas 
and dark matter present in this region. Considering the mass 
density profile given in Sect. |we can estimate the mass con- 
tent inside the effective radius of the Mrk86-C starburst pre- 
cursor. Assuming the most favourable case, i.e. i o =0° and 
a dark matter model with _Rdm=2 1,/2 x r, the density ex- 
pected at a galactocentric distance, r=1.0kpc, will be about 
9 atH cm -3 . Therefore, the mass content inside R c should be 



( |Fucntes-Massip 1997| ). Thus, we can conclude that the ve- 
locity dispersion in this region is sonic or slightly supersonic. 



7 SUMMARY AND CONCLUSIONS 

(i) The global velocity field in Mrk86 has a central an- 
gular velocity of about 34kms -1 kpc -1 with orientation 
PA~12° (rotation axis PA~ -78°). This velocity gradi- 
ent shows a progressive diminution towards outer galaxy 
regions, from 34kms -1 kpc -1 at the galactic centre to 
~10kms -1 kpc -1 at lkpc. The density profiles of the dif- 
ferent mass components indicate that the underlying stellar 
component dominates the total mass within its optical ra- 
dius. The velocity gradients measured for different distances 
along the galactic rotation axis indicate that the ionized gas 
is probably distributed in a inclined rotating disk, with in- 
clination in relation to the plane of the sky of about 50° . 

(ii) High velocity gradients of about 70kms -1 kpc -1 as- 
sociated with intense star forming regions have been ob- 
served. Both, a high mass concentration or a recent merger 
could be responsible for these local, steep velocity gradients. 

(iii) Our observations have revealed the existence of two 
bubbles, Mrk86-B and Mrk86-C. They present u cxp of 34 
and 17kms -1 and r bu bbie of 374 and 120pc, respectively. 
These structures are driven by supernovae and massive 
stars winds originated in two low metallicity (Z~0.1Zq) 
HII regions 11 and 8Myr old, with masses of 6.3 xlO 4 and 
5.3xlO 4 M . The H a emission is dominated by photoion- 
ization mechanisms, contributing at least with the 90 per 
cent of the total Ha luminosity. Moreover, the optical line 
ratios measured agree with photoionization mechanisms as 
the origin for their emission. 

(iv) In addition, we have studied the physical properties 
of the starb urst p recursor of the Mrk86-A bubble described 
by Martin (1998). The optical-nIR colours measured for 
these region are well described by a 14Myr old, low metal- 
licity (Z~O.25Z ), 1O 5 M massive starburst. 



t Turbulent velocity dispersion is denned as the line broadening 
obtained after subtracting instrumental and thermal broadening. 
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The predicted starburst ages agree quite well with the 
dynamical times measured if we assume a delay of ~4Myr 
between the starburst formation and the time when the bub- 
ble begins to inflate. In fact, the one tenth solar metallicity 
Leitherer & Heckman (1995) models predict an increment 
of about 1.5dex in the deposition rate of mechanical energy 
4Myr after the star formation burst. 

(v) The global velocity dispersion obtained for Mrk86-C 
does not trace the total mass of the burst. The turbulent ve- 
locity dispersion obtained for the Mrk86-C bubble is sonic 
or slightly supersonic, with <7t ur buicnt ~10 kms -1 . These val- 
ues do not show significant changes across the region. 
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